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Irrational constants in positronium decays 
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We establish irrational constants, that contribute to the positronium lifetime at 0{a) and 0{a^) order. In 
particular we show, that a new type of constants appear, which are not related to Euler-Zagier sums or multiple 
C values. 



1. Introduction 

Most of the multi-loop analytical calculations 
in quantum field theories have been dine for so- 
called single-scale problems. This means that the 
evaluated integrals are basicaly expressed as nu- 
merical constants up to a trivial scale factor. Ex- 
amples of such problems include almost all renor- 
malization group calculations, evaluations of the 
critical exponents, anomalous magnetic moments 
of the electron and the muon, matching calcula- 
tions in effective theories (e.g. HQFT, NRQFT) 
and many others. 

Usually analytical results involve the so-called 
Euler-Zagier (EZ) sums of the form 



E 



(±i)"i (±fc)"'= 



ni>n2>...>Tifc 



or more generally multiple polylogarithms 



E£]_ 



(1) 



(2) 



ni>n2>...>Tifc 



where zi, . . . Zk are some parameters and 
ai, . . . , flfc are positive integers. The sum ai -f 
a2 + ■ ■ ■ + cLk is called the weight in such a case. 



The above definitions include e.g. well-known 
irrationalities like C functions C,{a)^ Q{a^b), . . ., 
(poly)logarithms Lia(l/2), ln2, ... and "sixth 
root of unity" constants Ls^-'^'' (tt/S), Lsj'^'' (27r/3), 
.... There is no doubt, that by consideration of 
more complicated problems and in higher loops 
some mew constants will appear. As examples 
we can mention some elliptic integrals (see e.g. 

In this paper we concentrate on a very im- 
portant single-scale problem: the total width of 
positronium decay in QED. Positronium (Ps), the 
lightest known atom, provides an ultra-pure lab- 
oratory for high-precision tests of QED. In fact, 
thanks to the smallness of the electron mass m 
relative to typical hadronic mass scale, its the- 
oretical description is not plagued by strong in- 
teraction uncertainties and its properties, such as 
decay widths and energy levels can be calculated 
perturbatively in non-relativistic QED (NRQED) 
[4] with very high precision. 

Ps comes in two ground states, ^S{) parapositro- 
nium (p-Ps) and ^Si orthopositronium (o-Ps), 
which decay to two and three photons, respec- 
tively. 
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2. Orthopositronium 

In this section we are concerned with the hfe- 
time of o-Ps, which has been the subject of a 
vast number of theoretical and experimental in- 
vestigations. Its first precision measurement [5,, 
of 1968, had to wait nine years to be compared 
with first complete one-loop calculation [6] , which 
came two decades after the analogous calcula- 
tion for p-Ps [3 being considerably simpler ow- 
ing to the two-body final state. In the year 
1987, the Ann Arbor group [8 published a mea- 
surement that exceeded the theoretical prediction 
avalaible by ten experimental stantard deviations. 
This is so-called o-Ps lifetime puzzule triggered 
an avalanche of both experimental and theoreti- 
cal activities, which eventually resulted in what 
now appears to be the resolution of this puzzle. 
In fact, the 2003 measurements at Ann Arbor ^ 
and Tokio [10] 

r(Ann Arbor) = 7.0404(10) (8) ^is'^, 

r(Tokyo) = 7.0396(12)(11) fis-\ (3) 

agree mutually and with the present theoretical 
prediction, 

r(thcory) ^ 7.039979(11) ^ls-\ (4) 
The latter is evaluated from 
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(5) 



(6) 



is the LO result. The leading logarithmically 
enhanced 0{a^ In a) and 0{a'^ In^ a) terms were 
found in Refs. [12ll3j and Ref. [T3], respec- 
tively. The coefficients A = -10.286606(10) 
[6112115116117] , B = 45.06(26) [16], and C = 
—5.51702455(23) [18] are only available in numer- 
ical form so far. Comprehensive reviews of the 
present experimental and theoretical status of o- 
Ps may be found in Rcf. [T^ . 



Given the fundamental importance of Ps for 
atomic and particle physics, it is desirable to 
complete our knowledge of the QED prediction 
in Eq. ([5]). Since the theoretical uncertainty is 
presently dominated by the errors in the numer- 
ical evaluations of the coefficients A, B, and C, 
it is an urgent task to find them in analytical 
form, in terms of irrational numbers, which can 
be evaluated with arbitrary precision. In this Let- 
ter, this is achieved for A and C. The case of B is 
beyond the scope of presently available technol- 
ogy, since it involves two- loop five-point functions 
to be integrated over a three-body phase space. 
The quest for an analytic expression for A is a 
topic of old vintage: about 25 years ago, some of 
the simpler contributions to A, due to self-energy 
and outer and inner vertex corrections, were ob- 
tained analytically [21] , but further progress then 
soon came to a grinding halt. 




Figure 1. Feynman diagrams contributing to the 
total decay width of o-Ps at 0{a). Self-energy 
diagrams are not shown. Dashed and solid lines 
represent photons and electrons, respectively. 



The 0{a) contribution in Eq. dS]), Fi = 
ToAa/iT, is due to the Feynman diagrams where 
a virtual photon is attached in all possible ways 
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to the tree-level diagrams, with three real pho- 
tons linked to an open electron line with threshold 
kinematics. Such diagrams are shown in Fig. 1. 

After angular integration over three-photon 
phase space 



[dki][dk2][dk3]5{ki + k2 + ks - q) 



(7) 



we can rewrite the one-loop contribution to the 
width as (see [17] ) 



ma^ f dxi dx2 dx^ 
ri = TTTTj^ I — — d(2 - - a;2 - X3) 



SGtT^ J Xi X2 X3 



X [F{xi,X3) + perm.], 



(8) 



where Xi, with < a;^ < 1, is the energy of pho- 
ton i in the o-Ps rest frame normalized by its 
maximum value, the delta function ensures en- 
ergy conservation, and perm, stands for the other 
five permutations of xi,X2,X3. 

The function F includes dilogarithm and arc- 
tangent functions as given in [T7]. Ror illustra- 
tion, we just mention, that the above expression, 
after re-parametrization, consists of integrals of 
the following type 

1 

P{xi,X2,X3) f dy \n{xi + {1 - xi)y'^) 



Q{xi,X2,X3) J {I - Xi)x3, - Xi{l - X3)y'^ 



1 

P{xi,X2,X3) f dy \n{xi + {1 - xi)y'^) 



Q{X1,X2,X3) J X1X3 - (1 - Xl)(l - X3)y'^ ' 





with P, Q, P'Q' being some polynomials. 

The analytical integration of the above expres- 
sions is rather tedious and requires a number of 
tricks, e.g. expansion in series. Only a few inte- 
grals could be done strightforwardly, e.g.. with 
Mathematica or Maple. However, we established 
all irrational constants in terms of which the com- 
plete one-loop correction can be expressed. These 
include among others usual EZ sums up to weigth 
four, including e.g. 



In2. 



C(n) 



etc. 



and some additional constants of new type. At 
weight one, we have 

V2 — 1 

ln(i?) , where R — —= 

\/2 + l 

and up to weight four our basis includes the fol- 
lowing constants 



Li2 X 



Li 



1 
3 
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71 
Li4 (±i?) 



Li4 ( \ 



Lis (i?) 



Li4 ( -I 



Sl,2 (i?) , 
82,2 (±i?) , 



S1.3 (±i?) 

with Sa,b being the generalized polylogarithm 
1 

C_-| ^|a+&-l 

Sa.fc(a;)^ \ / — ln"-^Hn"(l-te 



l)\b\ 



-In'^-^Hn"^ 
t 



Unfortunately, not all integrals can be computed 
analytically. In more complicated cases, the inte- 
grations are not separated after expansion into in- 
finite series. We then rely on the PSLQ algorithm 
p6] , which allows one to reconstruct the represen- 
tation of a numerical result known to very high 
precision in terms of a linear combinations of a 
set of constants with rational coefficients, if that 
set is known beforehand. The experience gained 
with the explicit solution of the simpler integrals 
helps us to exhaust the relevent set. In order for 
PSLQ to work in our applications, the numeri- 
cal values of the integrals must be known up to 
typically 150 decimal figures. 

3. Parapositronium 

Let us now turn to the case of parapositronium. 
Its total width was recently measured to be [5S] 

rp(exp) = 7990.9Ais"^ (9) 

At present, the following radiative corrections 
within NRQED are available: 

me \ a / TT^ — 20 



— TT (-2Tr^ In a + A„) H In"^ a 



+ (^-Y + 101n2)lna 
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The first-order corrections were obtained in [7j, 
while the logarithmically enhanced terms were 
computed in [12ll3j . Here the constant Ap = 
5.12443(33) is known only numerically and 
our next goal is to establish the irrational con- 
stants that contribute to this quantity. 




in terms of special sums osf elliptic nature, 

n=l \2nJ > 

which we can call a^, and 02, respectively, 
and other sums 

f (-^6)" (1 11 (11) 

which we call 60 and h\. In ITUl) . stands for 
= 5-1(71 - 1) - 3S'i(2n - 1) + 251 (4n - 1) , 
with Sa{n) — X]j=i Vj° being a harmonic sum. 




Figure 2. Diagrams contributing to the decay 
width of p-Ps at 0(c?^ and their reduction to 
simpler master integrals. Dashed and solid lines 
represent massless and masive lines, respectively. 



This quantity recieves contributions from two- 
loop diagrams of e^e~ annihilation into two 
photons in threshold kinematics. However, the 
generic planar and non-planar diagrams (see 
Fig. 2, upper row) can be reduced via integration 
by parts to simpler integrals (Fig. 2, middle row). 
These, in turn, as we shall see, contain constants 
that are related to the sunset diagram (Fig. 2, 
bottom row) at very special kinematics, namely 
when the external momentum q is restricted by 

— —rr?. The sunset diagrams with such kine- 
matics have been considered in great detail in [T] . 
In particular the result for the sunset is expressed 



Figure 3. Example of vertex diagram J, con- 
tributing to the decay width of p-Ps. All lines 
have mass m. The dot on a line means the square 
of the propagator. 

Starting from pO|l and (fTT|) . one can construct 
sums of higher weights, e.g. 03, 0^2, ^3, etc. With 
such constructed sums, we evaluate more com- 
plicated diagrams, including vertexes and boxes. 
We illustrate it evaluating diagram J shown in 
Fig. 3. The resul is 

J=^C(3)-ia3-i«..-^63 (12) 

and a similar result follows for the box diagrams 
of Fig. 2. Formula (fT2|) shows the deep relation 
of the vertex diagram with the sunset diagram 
(in fact such relation follows from the differential 
equations). 

Concluding this section we want to mention 
that there are relations between the above sums 
and also their relation to the elliptic integrals has 
been found in [T]. 
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4. Conclusions 

Thus, we established the analytical structure 
of the results for the nex unknown corrections 
both for ortho- and parapositronium lifetimes. 
We found that new constants, that are not re- 
lated to the Euler-Zagier sums appear in both 
cases. 

We are grateful to G. S. Adkins for providing 
us with the computer code employed for the nu- 
merical analysis in [17] and M. Yu. Kalmykov for 
the fruitful discussions. 
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